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Abstract 

The spatially inhomogeneous large iV solutions to Kazakov-Migdal model are 
analyzed. The set of nonlinear differential equations is derived in the continuum 
limit. In one dimensional case these equations has a natural interpretation in terms 
of the dynamics of a Fermi gas. The multidimensional case seems to be inconsistent 
because of its instability related to the collapse of eigenvalues of the scalar field. 
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1 Introduction. 



A while ago, Kazakov and Migdal |l] proposed a lattice gauge model induced by a heavy 
scalar field in the adjoint representation of SU(N). The action of this model is the usual 
gauge invariant action without Yang-Mills term: 



S = -^N tr 



1 D 



u ($(x))-- ]T $(x)n v (x)$(x + v)nl 

1 v=-D 



X 



Although there is some problems with the induction of a physical QCD @, H, §J, |5| it is 
interesting to investigate the continuum limit of this model. In this paper we study a 
semiclassical dynamics of the density of eigenvalues of matrix : 

p(X,x) = ±tr6(X-$(x)). (1.2) 

An integral equations for p(X, x) can be obtained using the technique developed by Migdal 
H |5| . When the lattice spacing goes to zero a continuum limit can be constructed about 
the critical potential U cr ($) = D§ 2 . It is worth mentioning that the physical mass of 
scalar particles goes to zero in lattice units as the local limit is approaching in contrast 
to what is necessary for reproducing of QCD. 

A set of nonlinear differential equations is obtained in the continuum limit. In D = 1 
case these equations has a hydrodynamical interpretation and the translationally invari- 
ant solution is stable. When D > 1 the situation is qualitatively the other. The spectrum 
of excitations about the spatially homogeneous solution is always tachyonic. So the con- 
tinuum limit in D > 1 case seems to be physically unacceptable because of its instability. 



2 Saddle Point Equations. 

Due to gauge invariance all the matrices can be diagonalized by gauge transfor- 

mation. So, fixing the diagonal gauge and integrating over link variables one obtains the 
effective action depending upon - the eigenvalues of In the large N limit WKB 

approximation becomes exact. The semiclassical equation of motion reads as follows 

D 

-U^{X) + 2W{X,x) + J2 F u (\,x) = 0, (2.1) 

v=-D 

A varying along the support of p(X,x). The second term comes from gauge fixing deter- 
minant: 



W(\,x) = pjd£f^. (2.2) 



F u (X,x) is the logarithmic derivative of the Itzykson-Zuber integral: 

1 d 



F v (X,x) = lim ^^r-ln 1 2jn e JVtr *W M ^ nt 



(2.3) 
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Of course, it depends only upon the eigenvalue densities of and $(x + v). Using the 
Schwinger-Dyson equations for Itzykson-Zuber integral Migdal obtained the following 
dispersion relation determining F V (X, x) in terms of p(X,x) and p(X,x + v) ||: 

W(X, x + „) = / f arctan x _ F ^% (<x y (2-4) 

In spatially homogeneous case the equations ( |2.1| ), ( p.2|) and ( |2.4| ) were analyzed by 
Migdal in some detail and were solved exactly for the Gaussian potential by Gross . 
In the present paper we are interested in a spatial dynamics of the eigenvalue density. 



3 The Quadratic Potential. 

Before studying the general equations it is instructive to consider a more simple model 
with purely quadratic potential Uo($) = Ijlrn^ 2 . The saddle point equations can be 
simplified in this case due to the observation |7] that translationally invariant semi-circular 
distribution of eigenvalues solves (|2.1|), (|2.2p, (|2.4j). The semi-circular ansatz is useful in 
the case with a spatial fluctuations too: 



p(X,x) = Jp{x) -^p 2 {x)X 2 , F v =-f v (x)X. (3.1) 



Substituting ( |3.1|) into ( |2.2j) and ( pT4|) and doing the integrals we express W(X,x) and 
f u (x) in terms of /i(x). All the functions Uq(X), ^(A,^) and F v (X,x) are proportional to 
A with coeflitient of proportionality depending on the p(x) only, so the A-dependence in 
the equation ( |2.1| ) can be eliminated, that gives 



1 D 
p(x) = m 2 -- 

v=—D 



\ p(x + V) 



(3.2) 



This equation has translationally invariant solution J7| 



m 2 (D-l)±D^ml-4(2D-l) 
i± = • (3-3) 



It is interesting that ( |3.2| ) has no strongly fluctuating antiferromagnetic solutions. 

Now we are going to take the continuum limit of the equation Q3.2| ), so we rescale 
tUq = m 2 a 2 + 2D. At that moment a difference between D = 1 and D > 1 cases appears. 
Really, in the former case p + vanishes as the lattice spacing goes to zero. From ( p.3|) we 
see that p + scales as 2ma, so rescaling p(x) — > p(x)a and expanding (|3.2|) in a up to the 
second order we get for yjjj tr $ 2 (x)) = p^ 1 ^ 2 (x) = (j)(t), t = ix the equation 

0-m 2 + -L=O. (3.4) 
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This describes an oscillations with doubled frequency 2m about the static solution = 
(2m)" 1 / 2 . 

In multidimensional case /i_ vanishes in the local limit and scales as —m 2 a 2 /(D — 1), 
so after rescaling p,{x) — > p(x)a 2 we obtain from ( |3.2| ) the following equation (in the 
Minkowski space): 

□0-m 2 — = 0. (3.5) 

<? 

The effective potential for is unbounded from below, so the collapse of eigenvalues taking 
place in the naive (unregularized) continuum limit is inevitable. Translationally invariant 
solution = J—(D — l)/w? corresponds to a maximum of the effective potential, so it 
is unstable. 



4 The Continuum Limit for an Arbitrary Potential. 



Let us consider the one dimensional case first. The canonical scaling dimensions in ( |2.1|) , 
( ET2D and (fJJ) are recovered by the substitution A -> Aa" 1 / 2 , f -> ^a" 1 / 2 , U' Q (X) -> 
2Aa 1/2 + af7'(A), p(^,a?) -> p(^,x)a 1 / 2 . For F±(A,x) we can write F±(\,x) = AgT 1/2 + 
[v±(X,x) — W(X, x)\ a 1//2 + G±(A, x)a 3//2 . The first term in this expression is chosen so 
that ( f2.4| ) becomes the identity at the vanishing lattice spacing. The first and the second 
order terms in (12.41) reads 



c±(e,x) + 



(A-0 2 



^V(e^) 



A-e 



1 9 2 



2&r 2 



W(X,x). 



From (|2.1 ) we have 



v + (X,x) + ( A, x) = 
[/'(A) = G+(X,x) + G-{X,x). 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



The functions v + (£, x) and i>_(£, sc) differs only by sign, so we denote i>(£, t) = — iv+(^, 
s) = if_(^,x), where t = za;. After integration by parts in the l.h.s. of ( |4.1| ) we obtain 
the following equation for v 

dp 



9 , ^ 



0. 



(4.5) 



It is well known, that p(£,t) may be interpreted as the density of a Fermi gas in an 
external potential |§. So the equation ( f4.5|) has a natural interpretation as the hydrody- 
namical continuity condition, v being the velocity of a Fermi gas flow. The equation of 
motion for it is obtained as follows. First, differentiate ( |4.5|) in t, substitute the result in 
(|OD and integrate by parts in the r.h.s. Consequently integrating by parts the second 
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term in the square brakets. The last step is the elimination of G± by adding of the two 
equations ( |4.2| ) and using ( |4.4j ). The result reads 

r dp 



d, , 8 

dl {pv) + dt 



p[U + V 2 + -7rV 



u 



0. 



The present consideration is a generalization of a method used by Gross [ 
duce the well known fermionic solution || to one dimensional matrix model: 

1 



7T 



Po (0 = -J2E-2U(0. 



(4.6) 
to repro- 

(4.7) 



This is just the static solution to (|4.5p , ( |4.6| ). Constant E, the Fermi level, is determined 
by the normilization condition for po- 

In the acoustical approximation - p(£,£) = Po(£) + u(£,t), \u(£,t)\ <C Po(£), 



d 2 u 



d_ 



0. 



(4i 



The local velocity of sound, 7rpo(£), is equal to the Fermi momentum, as one might expect. 

Now let us turn to the multidimensional case. From the analysis of the Gaussian 
potential we learned that the density of eigenvalues scales as a 1 . So we write £ — ► £a _1 , 
A -> \a-\ p(£,x) -> p{£,x)a, U^X) -> 2D\ a - 1 +U'(X), F„(X,x) = Xa' 1 +v u (X,x) + 
[G U (X, x) — W(X, x)] a. It can be verified directly that this is the only way to obtain the 
sensible continuum limit of (2.1), (2.2) and (|2.4j ). All the steps in a derivation of the 
equations for p and v v are the same as in the one dimensional case. After the Wick 
rotation x D - 



ix°, Vd 



iv we get 
dp 



d 

+ (pv u ) = 



dx 



dx u dt; 

9 (pvn+p[U f + 2(D-l)W] + ^(pv' 



0. 



(4.9) 



(4.10) 



For the Gaussian potential the substitution of ( |3.1| ) with v v (£, x) linear in £ reduces these 
equations to (|3.5| ). 

Translationally invariant solution to ( |4.9| ), ( 4.10 ) is determined by the following con- 
dition 

Po(A) U'(£) 



dX 



(4.11) 



e-A 2{D-l) 

It coincides with the saddle point equation for one matrix model with a potential —U(£)/ (D- 
1). The generic solution to it is M, |9j 



MO = T JQ(& 



1 



(D-l) 



(4.12) 



where Q(0 is a polinomial of degree less than that of U'(£). However, this solution is 
unstable. The reason is that the r.h.s. of Q4.ll ) has a negative sign, so the eigenvalues of 
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$(x) are accumulated not in a minimum, but in a maximum of the potential. The insta- 
bility of the translationally invariant solution, of course, can be demonstrated explicitly 
by linearization of (fO|), ( |4.10| ) near po(0- The linearized equation reads 



n u + 2(D- 1] 



d_ 



Pop d\ 



u(\, x) 

7=X 



0. 



(4.13) 



Substitution u(£,x) = e " x «o(0 leads to the following eigenvalue problem 



Af{i) = -2{D-l)pj dX 



[Po(A)/(A)] / 



(4.14) 



There /(A) = l/p (A) /o A ^ M o(0- Operator A is self-adjoint with respect to the scalar 
product (/i,/2) = / <^A p (A)/i (A)/ 2 (A), all its eigenvalues being negative, because 

(f U Af 2 ) = -to - 1) f*dx Mi)fm - "° (A)/ y" mM^mm. 

J (4 - A) 2 



So all the spectrum of excitations is tachyonic. 



5 Conclusions. 

In multidimensional case we does not find such fine physical picture as in one dimension. 
Translationally invariant solution to the semiclassical equations of motion corresponds to 
a maximum of the potential, the fluctuations about it are unstable. Although we can not 
rule out the existance of a more complicated stable vacuum, we do not see the physical 
reasons for this possibility. So it seems that without the valuable changes, like an inclusion 
of fermions || |10j, Kazakov-Migdal model is unstable in the continuum limit. 
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